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ABSTRACT 

Nondestructive evaluation (NDE) method reliability can be determined by a statistical flaw detection study called 
probability of detection (POD) study. In many instances the NDE flaw detectability is given as a flaw size such as crack 
length. The flaw is either a crack or behaving like a crack in terms of affecting the structural integrity of the material. An 
alternate approach is to use a more complex flaw size parameter. The X-ray flaw size parameter, given here, takes into 
account many setup and geometric factors. The flaw size parameter relates to X-ray image contrast and is intended to 
have a monotonic correlation with the POD. Some factors such as set-up parameters including X-ray energy, exposure, 
detector sensitivity, and material type that are not accounted for in the flaw size parameter may be accounted for in the 
technique calibration and controlled to meet certain quality requirements. The proposed flaw size parameter and the 
computer application described here give an alternate approach to conduct the POD studies. Results of the POD study 
can be applied to reliably detect small flaws through better assessment of effect of interaction between various geometric 
parameters on the flaw detectability. Moreover, a contrast simulation algorithm for a simple part-source-detector 
geometry using calibration data is also provided for the POD estimation. 

Keywords: X-ray, flaw size parameter 


1. INTRODUCTION 

The paper provides modeling of X-ray flaw imaging contrast that is necessary to define the X-ray flaw size parameters. 
Commercial X-ray simulation software 1, 2 also model X-ray flaw detection and provide simulated images. The approach 
given here does not provide simulated X-ray images. Currently the commercial simulation software does not provide 
calculation of an X-ray flaw size parameter. However, similar to the approach given here, commercial X-ray simulation 
software can also be used to perform POD 3 studies. Comparison of results with commercial X-ray simulation software 
and use in POD studies is not covered in this paper. The model based X-ray parameter and simulated contrast are 
intended for use in POD studies and the approach can be referred to as a model assisted POD. The paper however does 
not provide any information on the POD analysis and limits to defining input quantities that are needed in the POD 
analysis process. 


2. NDE FLAW DETECTABILITY SIZE 

Crack location, crack orientation as well as crack opening influence the NDE flaw detectability. Flaw detectability is 
affected by choice of the nondestructive method, equipment capability e.g. sensitivity and resolution; and set-up 
parameters. Generally, a physical standard is used to calibrate and/or verify sensitivity and resolution requirements. Once 
performance of an NDE procedure in terms of the expected minimum results on a relevant physical standard is 
established, the procedure can be implemented for detection of real flaws. 

Flaw detectability of an NDE method can be given in terms of a reliably detectable flaw size. This is typically given as a 
flaw size with 90% probability of detection (with 95% confidence) and is denoted as a90/95 3 . A POD demonstration 
study is conducted where a set of flaws within the desired range of sizes is chosen for detection using a well documented 
and controlled NDE procedure. The a90/95 flaw size is calculated by analyzing the demonstration results using POD 
analysis software. It is important to remember that these estimations are dependent on number of flaws used in the 
demonstration. As the number of flaws in the POD demonstration increases, the a90/95 estimate tends to approach the 
true a90 value. It is customary to choose a definition of the flaw size parameter such that the POD increases with the 
flaw size parameter. Thus, the a90/95 flaw size is normally larger than the true a90 which is unknown. In some 
situations, especially where limited data is available, the a90/95 may not be assessed but instead a conservative value 
that is larger than the a90 is calculated. 


Page 1 


10/29/2013 



Surface crack length is the most often used flaw size parameter for assessment of a90/95 for dye penetrant testing and 
magnetic particle testing. Crack face area is the most often used flaw size parameter for assessment of a90/95 size for the 
ultrasonic testing. Crack depth-to-part-thickness ratio, a/t is commonly used for estimation of a90/95 in X-ray inspection 
of metals 4 . All flaw size parameters are based on physics of the NDE method. A monotonically increasing signal 
response is expected with increasing flaw size parameter value within a desired range of the parameter. 

3. X-RAY FLAW DETECTION ASSESSMENT APPROACH 

Many factors affect X-ray crack detectability. X-ray technique related factors are X-ray angle with respect to the crack 
plane, geometric unsharpness, detector plane angle, detector-part-source geometry, X-ray energy, exposure, detector 
type, detector sensitivity, detector noise, detector resolution, X-ray filter and single wall versus double wall part 
geometry. X-ray technique quality indicators (e.g. film density, 2% penetrameter hole detectability) are used to verify 
technique sensitivity and resolution. One of the factors is X-ray scatter which is dependent upon the energy of X-rays, 
part material and geometry. Moreover, the detector may have an effect on the X-ray scatter. The effect of X-ray scatter is 
not addressed in this work. The X-ray scatter and its effect is assumed to be low, uniform, and controlled. The scatter is 
considered to be one of the factors contributing to noise in the data. Crack, material and part related factors are: crack 
depth, crack opening, crack angle, crack location, flatness and roughness of crack faces, crack length; location of crack 
(e.g. in weld or parent material, surface roughness in the area of crack) and conditions that affect spatial noise (e.g. grain 
structure, scattering, part and weld geometry). 

Although, from fracture mechanics analysis point of view, it is desirable to have the NDE flaw size parameter to be only 
related to the crack major dimensions (e.g. length and depth), from NDE point of view, it is necessary to correlate flaw 
detection capability to both crack and key technique parameters. Currently, a90/95 crack depth-to-part thickness ratio 
(e.g. a/t) for film radiography in metal alloys is given as 70% for Standard NDE requirements 4 . The 70% number is 
arrived at by performing an X-ray flaw detectability study that meets the standard film X-ray requirements. 

We intend to build some models that relate selected factors such as the surface crack dimensions, orientation, part 
geometry, and technique parameters to film or image contrast for a simulated crack in a plate. Although the detector 
sensitivity, noise, X-ray exposure and X-ray attenuation coefficient affect flaw detectability, we assume that we could 
capture these factors in detector response calibration curves and as technique requirements. This allows us to focus on 
the geometric factors of flaw and source-part-detector geometry. A geometric flaw size parameter is intended to be 
monotonically related to the POD under these conditions. This approach assumes that the flaw size parameter and 
technique requirements including the part material requirements, the detector response calibration characteristics 
together provide a complete description of the X-ray flaw detection application and provide sufficient information for 
assessment of the NDE capability i.e. whether the flaw size under consideration can be reliably detected. A simulated 
contrast model provided here can do the same assessment if correlations are established between the simulated contrast 
and the actual contrast; and the simulation contrast and the POD. Applicability of these approaches would depend on 
validity of assumptions which affect accuracy of the models used. Therefore, experimental correlation with the signal 
response or POD is necessary to validate the models. 

4. OUTLINE OF THE MODEL DEVELOPMENT 

The paper is fairly lengthy in terms number of equations. Therefore, some explanation of steps taken to build the X-ray 
parameter and contrast model are necessary. Many quantities and some modeling approaches must be defined before the 
X-ray parameter and contrast can be defined. Initially, basic definitions and descriptions of the crack shape, geometry of 
the set-up; X-ray attenuation, X-ray intensity, exposure, film density, image intensity, and contrast are given. The 
definitions are followed by normalized detector response from a step wedge. The detector can be a film or a digital 
detector. Linear and non-linear detector responses are provided as a function of ratio (a/t) of a step wedge. The models 
include a noise term. Definitions of the contrast and normalized contrast are given. These definitions and models are 
necessary for development of the X-ray parameter in later sections. Definitions of relative noise and contrast-to-noise 
ratio are given as these quantities are useful in the POD study. In order to model a detector, a definition and model of the 
modulation transfer function are given. The detector contrast and normalized contrast that account for effect of the 
modulation transfer function are provided as a function of the X-ray parameter. Next, three models of the X-ray 
parameter are developed. A computer simulation is provided to illustrate the values of the X-ray parameter and contrast. 
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Finally, a computer application is developed to calculate the X-ray parameters and contrasts. Explanation of use of these 
quantities in the X-ray NDE POD study is provided. The paper, however, does not provide any experimental data in the 
use of these quantities in a POD study. 

Fig. 1 provides four separate blocks of the model development and general progression of the model development 
between these blocks. The first block lists basic quantities that have been defined. These quantities are necessary in the 
development of model. The second block provides a flow chart for modeling of the X-ray parameter. The third block 
provides a flow chart for modeling of the detector response. The fourth block provides a flow chart of modeling of MTF 
and detector response accounted X-ray parameter and contrast. 




Fig. 1 : Flow chart of the model development 

5. MODELING CRACK SHAPE AND SETUP GEOMETRY 

Here, we take a single open-to-surface rectangular slot as the crack in a plate. The slot depth is denoted by a and width is 
denoted by W. The plate thickness is denoted by t. The slot is normal to the plate surface. The imaging detector can be a 
film, computed radiography (CR) screen or any of the variety of digital detector panels or cameras. See Fig. 2 for profile 
of the slot. The X-ray beam is in the plane of figure. The slot length and plane of detector plate are normal to the plane of 
figure. We assume that an estimate of crack width is available. The material thickness under the slot is denoted by b. The 
X-ray flaw size or X-ray parameter models will be developed for this part- source-detector geometry. 
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Parallel X-rays 



Fig. 2: Part, slot, detector geometry and selected X-rays and X-ray shadow profile 

We also use a step wedge 5 made from the plate material for the signal response calibration. See Fig. 3. Thickness of the 
thickest step of the wedge is same as the plate thickness. Thickness of each step is then denoted by b and depth of each 
step from top of the highest step is denoted by a. 

We also use line-pair resolution targets 5 to characterize the modulation transfer function of the detector as applicable to 
the X-ray test set-up. See Fig. 3. The plate with a slot, the step wedge and the line-pair resolution target are subjected to 
the same X-ray inspection set-up. Slot depth a, material thickness under slot b and plate thickness t are related by, 

b = t - a or a = t - b. (1) 

Since the flaw detection is based on image contrast of the flaw indication, we intend to formulate the flaw size parameter 
that relates to the image contrast. 



Fig. 3: An aluminum step wedge and a line pair resolution target made from 0.05 mm thick lead shim 


6. X-RAY ATTENUATION AND FILM DENSITY 

Let us first consider a radiography setup where a plate of thickness t is radiographed. We assume that the X-rays are 
normal to the top surface of the plate and the detector is in intimate contact with the bottom of plate. If intensity of X- 
rays entering the plate is I 0 and intensity of X-rays exiting the plate is I t , then the two intensities are related by 6 7 ’ 8 , 

(1) I = I,e~ M o \nl,=\nl 0 -at (2) 

where, 

Jo = X-ray intensity or exposure incident on top side the of plate, 
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I t = X-ray intensity or exposure exiting from bottom side the of plate, 
a = X-ray attenuation coefficient. 

The above equation is valid for monochromatic X-radiation and is assumed to be approximately valid for X-ray tubes 
used in nondestructive evaluation. X-ray intensity is measured in rem per unit area per unit time. Exposure or absorbed 
energy dose is measured in rem 6 ’ 7 . However, in industrial radiography, the exposure is defined differently. It is defined as 
the tube current times the shot duration 6 . The source energy, source-to-part and part-to-detector distances shall be noted 
as change in the source-to-detector distance would affect the exposure. We use the X-ray intensity and exposure (tube 
current times the shot time) interchangeably because the X-ray intensity multiplied by area and time is proportional to 
the exposure in a given set-up 6 . 


X-ray film density is measured by noting fraction of light transmitted through the film. The film density is defined by the 
following expression 6 , 


d = log 10 ^, 
A 


(3) 


where, 

if = intensity of light incident on the film, 


if = intensity of light transmitted through the film, 
d = film density. 


The film density correlates better with logarithm of the X-ray exposure. 

7. IMAGE CONTRAST 

The result of an X-ray inspection is a digital image or film with an image of the slot. The slot is identified in the image 
due to image contrast. Contrast or net detector response is defined as 6,7 , 


where, 

J a = detector response registered at the slot image location, 

Jo = detector response registered in the surrounding region of plate thickness t (or at a = 0). 

In the above equation we assume that the slot indication has a higher numerical value for image intensity or film density 
than in the surrounding region. See Fig. 4 for an example of the net detector response for a film obtained on a step wedge 
as a function of step depth-to-part thickness ratio (a/t). Here, the data-points have been joined by a smooth curve. 
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a/t for Step Wedge 
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Fig. 4: Detector response function for film 


Here, we introduce normalized contrast which is based on the net detector response. Net detector response due to a slot 
or step of depth nn a plate of thickness t is given by, 


S , - J , J o, (5) 

where, 

J t - detector response at a - t (i.e. plate is absent). 

By definition, the net detector response for a = 0 is zero. S 0 = J 0 — J 0 = 0. The normalized contrast or normalized 
detector response, as defined here, is given by, 



• ( 6 ) 
See Fig. 5 for an example of the normalized detector response or normalized contrast for a digital detector obtained on a 
step wedge as a function of step depth-to-part thickness ratio (a/t). It is assumed that an adequate exposure is used so that 
the signal-to-noise (S/N) ratio is high. 



Fig. 5: Normalized detector response function for a digital detector 

8. LINEAR DETECTOR RESPONSE MODEL 


Let us consider a linear model of the detector response first. We assume that the detector contrast is approximately linear 
with logarithm of the exposure 6,7 . Using eq. (2) for log of exposure, the detector response is modeled here as, 

J 0 =q 1 (\nI 0 -at-q 3 ) + q 2 +S , (7) 

where, 

Jo = detector response due to X-ray exposure on a plate of thickness t , 

C[i = detector sensitivity, 
q 2 , q 3 = detector coefficients, 

3 = detector noise. 

Detector noise is assumed to be normal with zero mean. Standard deviation of the detector noise is denoted by (J § . We 

assume that the detector response is averaged in the background. If an average value is used for J 0 , then the noise term 3 
is not used. The detector does not respond until it receives a minimum level of exposure. Also, the response is saturated 
above a certain level of exposure. Thus, the linear model is valid only between these extremities. The digital detectors 
have a higher dynamic range and respond to relatively low exposure as well as high exposure with a linear relationship 
throughout the range. The X-ray film has a lower dynamic range compared to the dynamic range of an amorphous - 
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silicon digital detector. Relationship between the film density and exposure is non-linear. Relationship between the 
digital detector response and exposure is approximately linear. For convenience we would omit the noise term in the 
following five equations. By substituting step thickness b and t separately in eq. (7), various detector response 
expressions are derived as, 

J„ =q l (\nl 0 -ab-q 3 ) + q 2 
J a = q x (aa) + J 0 , and 


( 8 ) 

(9) 


J t q l (^at^j + J 0 

Here, we assume that the measurement is done on an indication of a step in a step wedge or on a wide slot in a plate. By 
rearranging eq. (9), the contrast or net detector response at a step is given by, 

S a= J a~ J Q =q,CCa 

By rearranging eq. (10), the net detector response for the plate is given by, 

S t = J t — J 0 = q x at 

By substituting eq. (11) and (12) in eq. (6) and neglecting noise, the ideal normalized contrast is given by, 


( 11 ) 

( 12 ) 


C = -^- = a/t . (13) 

However, in practice the contrast measurement would be affected by noise. Following expression indicates effect of the 
detector noise 8 (fi, a) in measuring contrast. 


q x aa + 8 
q x at 


(14) 


Here, the noise term in the denominator of the above expression is neglected by assuming averaging of intensities over 
appropriate areas in measuring S t9 


Contrast and normalized contrast are assumed to relate to flaw detection ability. The above equation is independent of J 0 
but assumes that logarithm of the exposure provides a response that is within the linear response range. The above 
equation implies that higher step (or slot) depth-to-part-thickness ratio provides higher normalized contrast. 


Contrast-to-noise ratio 9 is given by, 




(15) 


J S ^8 

The above equation implies that higher material attenuation coefficient, higher detector sensitivity and slot depth provide 
better signal-to-noise ratio. A high contrast-to-noise ratio (e.g. > 3) is desired. Relative noise, as defined here, is given 
by, 


R noise = y- • ( 16 ) 

*'0 

Detector noise is used in noise analysis which is used to supplement the POD analysis. Detector exposure factors and 
contrast-to-noise ratio are better addressed in the contrast than in the normalized contrast. Normalized contrast is related 
to flaw detection and it is less dependent upon exposure, X-ray attenuation coefficient and detector characteristics. Later, 
ratio a/t is extended for narrow width slots by defining the X-ray parameter P that is related to the flaw contrast. 


9. NON-LINEAR DETECTOR RESPONSE MODEL 


An example of a non-linear model is given in Fig. 4. By neglecting detector noise, and assuming a typical nonlinear 
detector response as shown in Fig. 4, a general non-linear model can be given as, 

4 =/„(«/» Wo 


•W.W+A 
/„( 0)=0 


(17) 

(18) 

(19) 

( 20 ) 
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where, 

J a = detector response due to X-ray exposure to plate of thickness b , 
fn( a / t ) = a function or dependent variable of a/t, 

J 0 = detector response when a - 0 in a part thickness of t. 


The ideal normalized contrast is given by, 

c= [M t ) 

' m ' 


( 21 ) 


However, in practice the contrast measurement would have the effect of noise. Following expression indicates how 
contrast is affected by the detector noise. 


f{a/t) + S_ 

m 


( 22 ) 


9. MODULATION TRANSFER FUNCTION 

Modulation transfer function (MTF) as a function of line-pair per frequency (or density of an image target is commonly 
used to characterize resolution of an imaging array or film type detector. A line-pair in an imaging target is defined by a 
dark line and a light line of equal width L placed next to each other 6,7 . The line-pair width in the target is given by, 

Li p = 2 L (23) 

Line-pair frequency measured in line -pair per mm is denoted by w tp . The line-pair frequency and the line-pair width are 
related by, 

w i P = y- • (24) 

L i P 

We assume that the width L ip is given in millimeter. A 2D digital detector array panel has clearly defined pitch between 
pixels. The pitch or the detector pixel width can be taken as the detector width L d . 


The line-pair resolution is also applicable for X-ray inspection. In the X-ray inspection, a line-pair target is made from a 
thin lead shim. See Fig. 3. Many pairs of shim-gap with same width are grouped so that the image width of any line-pair 
group is longer than the detector pixel size. 

Modulation transfer function can be calculated by imaging the line pair resolution target and using the following 
equation. 

M = 2) — A (25) 

J t +J 0 

where, 

Jq= minimum film density or image intensity of the shim indication, 

][- maximum film density or image intensity of the gap indication. 


In the above equation we assume that the gap indication has higher numerical value of image intensity than the shim 
indication. MTF is plotted and modeled as a dependent variable of w tp or L ip of the imaging target at the detector. When 
an indication is very wide i.e. more than ten times of the detector element size, the registered intensity is altered very 
slightly due to slight change in the indication width. The maximum MTF is given by, 


Normalized MTF is given by 10 , 


M =2^0 

max j j 

J t + J 0 


M. = 


M 


M 


(26) 


(27) 
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( 28 ) 


Sum of the modulated intensities is same as that of unmodulated intensities. Therefore, 

3 1 A) — 3 1 • 

Therefore, by combining eq. (26) and (27) the modulated intensities relate to the unmodulated intensities as, 
Sum of the modulated intensities is same as that of unmodulated intensities. Therefore, 

J a + ^0 = J a + ^0 • 

By combining eq. (27) and (29) the modulated intensities relate to the unmodulated intensities as, 


(29) 

(30) 

(31) 


The modulated signal response at the slot for a line-pair resolution target is given by, 

/,+r 0 + M,(r,-r 0 ) 

“ 2 

The above equation is true for imaging a line-pair target. MTF may also be a function of angle of the incident radiation. 
If MTF is used in the analysis, then detector unsharpness U d should not be used. The normalized MTF is modeled as a 
function g n with line -pair as the dependent variable. 

M. =g,( w i,) or ( 33 > 

M„=g,(l/L, f ). (34) 

A curve may be modeled in the MTF plot. For example, the modulation transfer function may be modeled as, 

SnM = a l e a2W,p +a/ iW \ (35) 

where, a 2 , a 3 and a 4 are the coefficients. 


Examples of assumed MTF for two different detectors are given in Fig. 6 and 7. Here, we choose a line-pair resolution 
corresponding to an MTF value of 0.4 as the line -pair resolution (L d ) of the detector. The detector resolution is assumed 
to relate to the pixel size of a digital radiographic detector or average grain size of a radiographic film. Another MTF 
value (e.g. 0.2) can be chosen. The chosen MTF value allows us to quantify the resolution of a detector in terms of line- 
pair per frequency. 

MTF value for a film is 0.4 for a 20 lp/mm target as shown in Fig. 6. 
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Fig. 6: MTF used for film in the simulation 


MTF value for a digital detector is 0.4 for a 4.2 lp/mm target as shown in Fig. 7. 
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Detector Element Size: 120 micron Detector Element Size: 120 micron 
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Fig. 7: MTF used for a digital detector in the simulation 


By substituting the flaw size parameter P, defined later, in place of a/t in eq. (21) and multiplying by the normalized 
MTF we get MTF accounted normalized contrast parameter. 

c =M = 2 (, \fjfl ,,,, 

If the detector response is not known, then the normalized contrast parameter based on the X-ray parameter and MTF is 
given as, 

C„,=M n P = g n (yi 1 p )P. (37) 

Here, we rewrite the above equation as, 

p K s=M.p=g,(yh)p- < 3 *> 

The above quantity can also be considered to be modulation accounted X-ray parameter. The above two equations relate 
to the normalized contrast of a single slot when modulation effect of the detector is accounted. Modulation accounted 
contrast is given by, 

s MAP] =Mj n {p)=g n (\m P )f n {p) (39) 

Similarly, it is assumed that the shadow contrast f n (P) is uniform within the shadow width. 

11. MODEL BASED ON EQUIVALENT X-RAY PARAMETER 

Visual perception of the slot indication is primarily influenced by the image contrast. The perception is affected by 
background intensity to a lesser extent. Visual perception is not good if images are too dark or too bright. Therefore, in 
film evaluation, the viewing light intensity is controlled by a backlight illuminating the film. In digital images, the image 
processing software allows changing contrast (and brightness) by rescaling within a chosen viewing intensity range. Base 
intensity should also have good signal-to-noise ratio. Typically, contrast-to-noise ratio of greater than 3 is needed to 
detect flaws reliably. 

Ratio a/t is a measure of the ideal normalized contrast. Therefore, ratio a/t can be considered to be a fraction of the X-ray 
pathlength passing through void of the slot. Ratio a/t can be described as the “void pathlength ratio”. During inspection, 
an X-ray penetrameter made from the plate material, is placed on top of the plate. A 2-2T hole of the strip 
penetrameter 10 provides a void pathlength ratio of 0.0196 or approximately 2%. 

If the 2-2T hole or the penetrameter edge can be identified in an X-ray image, it does not imply that a crack with 2% 
depth can be imaged or detected in the radiographic image. This is due to combined effect of the crack related factors 
and technique related factors. Some of the technique related factors include incident angle of X-rays, geometric 
unsharpness and resolution of the recording medium. Some of the crack related factors include crack depth, length, 
opening, straightness of crack faces, and roughness of crack faces. 
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The model uses a parameter based on the spatial variation of the X-ray pathlength ratio through the void as mapped on 
the recording medium. The contrast is also affected by geometric unsharpness and other factors. Therefore, the X-ray 
parameter is based on assessing average effect of the X-ray pathlength ratio that also includes the effect of geometric 
unsharpness in the image of slot. Spatial variation of the pathlength ratio through the slot in combination with 
unsharpness effects is called the X-ray shadow profile here. We do not intend to provide an algorithm to compute 
shadow intensity for every point in the shadow profile. We assess an equivalent intensity for the slot shadow in the X-ray 
parameter. The maximum possible value of the X-ray parameter is a/t. 

Geometric unsharpness causes edges of the shadow to be blurred. Zone of the blurred edge of the shadow is called 
penumbra. Middle (or inner) portion of the shadow is called umbra. Width measured between the outer edges of 
penumbra provides the outer shadow width. Equivalent shadow width (L e ) is estimated as a centroid based or simple 
average of the outer and inner shadow widths. Two times the equivalent width of shadow is used as the line-pair width in 
computing the MTF. We use the X-ray parameter and associated equivalent width in evaluation of the normalized 
contrast and contrast. 

In reality, X-rays may have an oblique angle to the crack face. The crack opening may be very narrow. These factors 
would decrease the X-ray parameter below its maximum possible value of a/t. The X-ray parameter is used in the MTF 
accounted normalized contrast and contrast as follows. 

Q. /lPj = .?„(l/2A )y^y and (40) 

S„ m = 8n(yu*)fn(P)- ( 41 ) 

Consider a special case of an ideal detector where the detector response is linear, 

f n {P) = kP and /„ (l) = k . (42) 

Here k is the detector sensitivity constant in a given inspection set-up. Note that k also depends upon the X-ray 
attenuation coefficient. Substituting eq. (42) in eq. (40), we get, 

C M ,f {P) =g n (VZLe)P- (43) 

We would rewrite the above equation as, 

Pm,/(p) = S n (V2L e )P • (44) 

The above quantity can also be considered to be the modulation accounted X-ray parameter. The above equation relates 
to the contrast of the line -pair resolution target but also relates to normalized contrast due to a single slot. 

S M , m = 8n{V2Le)kP. (45) 

Contrast for a single slot is affected by the detector pixel size. The above two equations indicate that, when the net 
detector response is linear with a/t , the contrast and the normalized contrast are proportional to a product of the X-ray 
parameter P and MTF value for the shadow width. A high MTF value is obtained if the detector has high resolution in 
comparison to the shadow width L e . A high value of detector sensitivity k , near the highest possible value for the 
technique, is desired to maximize the contrast. Normally, radiographic film sensitivity is high for a film density of 2.5 to 
4. The radiographic technique (energy, filters, duration, density) can be controlled such that a desired value of the 
detector sensitivity k can be obtained. 

The above equations indicate that the X-ray parameter relates to the contrast which is assumed to relate to the flaw 
detection. Therefore, the X-ray parameter relates to the flaw detectability or POD, provided the detector sensitivity and 
detector resolution are controlled or standardized. The technique sensitivity can be assessed by radiographing and 
evaluating a resolution target and step wedge. 

We would further develop the X-ray parameter so that it relates to some more set-up geometry related parameters, 
including source-to-part distance, part-to-detector distance, X-ray angle, source size, slot width, slot depth, geometric 
and other unsharpness. 

12. MODEL LEVEL 1: PARALLEL X-RAYS 
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In developing a model for the X-ray parameter, we start with the simplest configuration. Here we assume that the source 
provides parallel rays in level 1 model. A crack is modeled as a slot described earlier. See Fig. 2. Two modes (e.g. mode 
1 and mode 2) are defined here depending on interaction of a parallel X-ray beam with cross sectional geometry of the 
slot. 


Mode 1 is defined by the following condition, 

W>atan/?, (46) 

where, 

/? = angle of X-rays with respect to the normal to the part surface. 


When this condition is satisfied, some X-rays entering the slot opening on top side pass through flat bottom of the slot. In 
mode 1, P max ,i is given by, 

P m ^=a/t, (47) 

where, 

P max, i = peak value of the X-ray shadow profile. See Fig. 2 for a schematic of the X-ray shadow profile. 

Subscript “1” is used to indicate the model level. A slot of width W casts an image with length given by L 2 i, 

Lu=W + atanj3 (48) 

L 2 i is also called outer width of the shadow. Slot indication contrast is high in the center (zone called umbra) and it 
tapers to the side (zone called penumbra). Width of the center or inner portion with higher contrast is given by L u . 

Z 1I = W -a tan (3 . (49) 

L lsl . is also called inner width of the shadow. Mode 2 is defined by the following condition. 

W<atan/?. (50) 

When this condition is satisfied, no X-rays entering from top of the slot opening pass through flat bottom of the slot 
opening. The maximum X-ray shadow intensity is given by, 

D w 

P m ax,l=- (51) 


t tan p 

In mode 2, outer width of the shadow is given by, 

L 2l =W + a tan /? 

In mode 2, inner width or width of the center portion of shadow is given by, 

L, a = a tan fi- W 


(52) 

(53) 


There is a smooth transition between mode 1 and mode 2 as the X-ray angle increases. This is evident by a common 
point shared by both modes at, 

W = a tan/?. (54) 

Here, we generically define the two shadow widths as, 

L hn = umbra or inner shadow width. Subscript “w” indicates the model level number 
L 2j n = total or outer shadow width. Contains both umbra and penumbra shadow zones. 


For both modes, the inner shadow width is given by, 

L, 4 =|atan /? —W\ . 


(55) 


Human detection of a fine crack is assumed to be dependent not only on the contrast associated with peak of the shadow 
but also on indication widths L x and L 2 . This assumption will not hold if L x and L 2 are about same or value of X-ray 
parameter is too high (e.g. possibly greater than 10%). But barring these extremities, the area under the shadow profile 
can be used to assess equivalent indication intensity. The X-ray parameter is given by, 


P = P 


r 2 

A +L 2 


(56) 


Subscripts indicating the model number have been dropped to indicate the general formula. Subscript “c” is used to 
indicate that the X-ray parameter is based on the centroidal weighing or average i.e. centroidal distance (height) of 
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shadow profile in order to assign a single value related to contrast perception. See Fig. 2. P c is defined as two times the 
height of centroid of the shadow profile in Fig. 2. The X-ray parameter is not used to calculate an estimated image 
contrast but is used to calculate an “equivalent” contrast that may be related to the actual contrast. Here, we can also use 
an approach to divide the area under the contrast curve by the average width of the indication to obtain a quantity related 
to contrast. While such approach may also work, the centroid is believed to provide a better weighing of the area under 
the contrast curve. 


P c is equal to twice the centroidal height of the indication. Thus, if L x and L 2 are equal then P c = P max and if the shadow 
profile has a triangular shape (L x = 0), then P c = (2/3 )P max . By substituting eq. (47) in eq. (56), for mode 1 we get, 


PcX 



r 2L, a + L l l ^ 
v A.1 + k,l J 


(57) 


By substituting eq. (51) in eq. (56), for mode 2 we get, 


p ->-\ 


f 


w 


t tan p 


2 A,1 + k,l 

A.1 k,\ 


(58) 


13. MODEL LEVEL 2: FINITE SIZE SOURCE AND UNSHARPNESS 


Here we begin work on a slightly more sophisticated level 2 model. In this model, we account for the unsharpness in a 
simple way. We define unsharpness as follows. 

U n =U g +U a (59) 

Here subscript “g” is used to indicate the geometric unsharpness. Subscript “o” is used to indicate unsharpness due to 
other effects including film (subscript “/’), forward scatter (subscript “s”) and cumulative of other miscellaneous effects 
(subscript “m”) 

U 0 =U f +U,+U m (60) 


Although, unsharpness components are identified, in most cases only geometric unsharpness is known. Other factors 
contributing to the unsharpness can be considered to be part of the conditions that would be controlled. Under these 
assumptions for mode 1 and 2, the inner shadow width is given by, 


ki = 


atmfi-W\ 


^ + (^2 ^ 


d j 


-U, 


(61) 


where, 

dj = source to part distance and 
d 2 = part to detector distance. 


The above equation is derived by applying geometric enlargement to the inner shadow width of model 1 and then 
subtracting the net unsharpness. The outer shadow width is given by, 


L 22 =(atanf3 + W) 


^ d 2 ^ 


+ [/ 


(62) 


The above equation is derived by applying geometric enlargement to the outer shadow width of model 1 and then adding 
the net unsharpness. 


The average indication length for model 2 is given by, 

A , 2 + ^ 2,2 


L e,2 = 


(63) 


X-ray parameter for mode 1 in model 2 is given by, 
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2 


22,.,+i,,,) 

Aj + O 

T 

v t y 

v A, 2 + A, 2 y 

v A, 2 + A, 2 y 


( 64 ) 


The above equation is similar to corresponding equation for model 1 , except the increase in length of indication due to 
geometric enlargement and unsharpness is accounted by proportionately reducing the X-ray parameter. Similarly, X-ray 
parameter for mode 2 in model 2 is given by, 


2 

r w j 

2A, 2 + A. 2 ) 

A,, + 0 

3 

v /tan/? y 

v A, 2 + A, 2 y 

v A, 2 + A, 2 y 


(65) 


14. MODEL LEVEL 3: FINITE SIZE SOURCE AND UNSHARPNESS 

Here we begin work on model 3 which will be slightly more sophisticated than model 2 but provides the same quantities 
with slightly different equations. Similar to model 2, model 3 is derived from model 1. In model 3, the unsharpness is 
handled differently. Unsharpness is handled as mainly contributed by the geometric unsharpness. Geometric unsharpness 
is given by, 

Sd 2 

U=—^, (66) 


where, 

S = width of the X-ray source. 

Here, we make an assumption that the slot can be replaced by an equivalent shadow in the shape of a rectangle with 
height given by, 


P =P =-P 

1 x, 3 c, 1 1 max,l 


2 i-Lu+L,,) 


3 — (2,,,+Z,,) ' 


and the equivalent width is given by L x . 


J _ 3 (a,! + ^2,l) 

" 4 (2Z,J + 


(67) 


( 68 ) 


L x is derived by keeping area under the shadow profile same as in model 1. The above equation provides a length 
inbetween L u and L l x The length is approximately given by L u + (2/3 )(L 2 ,i - L u ). 

The simulated slot shadow with width W = L x3 is then projected onto the imaging plane using X-rays from the source. 

In this situation, the outer and inner widths of the shadow are, 


A, 3 


k 3 = 


At, 3 (A + <A) 
d x 

A, 3 (A + A ) 


d \ 


+ U ' and 


-U 


(69) 


(70) 


We assume that area under the shadow profile is same between a parallel beam source and a source with known size. 
Therefore, 


Aiax,3 ( Al,3 + A, 3 ) — Aiax,l ( A,1 + A,1 ) ’ 0r 

A,1 + A,1 


P = P 

max, 3 max,l 


A, 3 + A, 3 


Using eq. (55), X-ray parameter is given by, 

p --P ^ A,3 + A, 3 

c ,3 ^ max, 3 t- 

* Al,3 A, 3 


(71) 

(72) 

(73) 
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Equivalent shadow width is given by, 


L = 


aW 1 


t P 


or 


c,3 


3 (A, 3 + ^2,3) 

^ A, 3 + ^2,3) 


(74) 


(75) 


The equations for model 3 are slightly different for equivalent width than those of model 2. One of the two models (e.g. 
model 2 or 3) is likely to be slightly better depending upon magnitudes of geometric unsharpness, other unsharpness and 
geometric enlargement. 

15. DETECTOR ORIENTATION AND OTHER FACTORS 

Here we consider an angle a between the detector plane and plane of the plate. Angle between the detector and the plate 
increases width of the shadow on the detector. The corresponding shadow lengths are. 

L va =L v 3 sec a and (76) 

Z^ 4 =7^3 sec a. (77) 

The larger shadow width is advantages in detecting smaller flaws but the relative noise increases. Moreover, the detector 
may be less sensitive to obliquely incident X-rays. The X-ray parameter is unaffected. 


The average shadow width is given by, 


p = p 

1 c, 4 1 c , 3 • 


L eA= L e,3 SQCa - 


(78) 


(79) 


Equivalent width is useful to model detector resolution effects. We also consider the unsharpness U d at the detector 
screen, especially if X-ray radiation is not normal to the detector. Detector unsharpness may be a function of angle of the 
incident radiation. Equivalent width of the shadow at the detector and the X-ray parameter are given by, 

L e,5= L eA+U d and (80) 


P = P 

1 c, 5 1 c, 4 




~^e,5 


( 81 ) 


16. SIMULATION RESULTS 

A set of rectangular slots are taken for simulation. Two detectors with corresponding discrete modulation transfer 
function g(wj p ) and the signal response function f(a/t) are used. One detector is assumed to simulate a film and the other 
is assumed to simulate a digital detector. It is assumed that the detector responses are obtained by radiographing a step 
wedge with maximum step thickness of t. The detector response is plotted as a function of a/t. 

Fig. 4 shows the discrete net detector response function ( S a ) for a film used in the simulation. The function is given in 
film density units. The base density J 0 is 0.45. Fig. 5 shows a detector response function as normalized contrast (C) for a 
digital detector. The digital detector response is given as the normalized response between 0 and 1 and is assumed to be 
equal to the a/t of the step wedge used. 

X-ray parameter P c 3 values for various combinations of slot width and depth for two different angles of incidence are 
calculated. Three part thicknesses are chosen. These are 0.28, 1.128 and 2.4 mm. The X-ray parameter values are plotted 
as a surface. Three surfaces, one for each part thickness, are shown. Fig. 8 shows the X-ray parameter values for a 
normal angle of incidence. The source size is chosen to be 1.2 mm and the source to part distance of 60 cm and part to 
detector distance of 10 cm. 
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X-ray Angle = 0 degree, Source Size = 1 .200 mm, Detector Size = 0.200 mm 



Fig. 8: X-ray parameter, source =1.2 mm 


The above figure illustrates that the X-ray parameter increases with increase in crack width or depth and decreases with 
increase in the plate thickness. 

Effect of the oblique X-ray angle is illustrated in Fig. 9. 


X-ray Angle = 5 degree, Source Size = 1 .200 mm, Detector Size = 0.200 mm 



Slot Width, mm 0 Slot □gp^ mrT1 


Fig. 9: X-ray parameter, source =1.2 mm 

The above figure illustrates that the X-ray parameter decreases with increase in angle of incidence of X-rays. Effect of 
the source size is illustrated in Fig. 10. 
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Fig. 10: X-ray parameter, source = 0.2 mm 

Smaller source size is used in Fig. 10 compared to that used in Fig. 8. It illustrates that the X-ray parameter increases 
with decreasing source size. Many other interesting combinations of input parameters and resulting charts are possible 
but are not presented here. MTF accounted X-ray parameter is shown below. 

X-ray Angle = 0 degree, Source Size = 0.200 mm, Detector Size = 0.125 mm 


5 
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Fig. 11: X-ray parameter, source = 0.2 mm, detector resolution 4 lp/mm 



This parameter P C} 3 >M accounts for many source, part and detector geometry factors. It does not account for internal 
scatter, and detector response. The detector response is dependent upon factors such as material attenuation coefficient, 
X-ray energy and exposure. Typically, penetrameter sensitivity of 2% is demonstrated. This will be equivalent to P c ,3,m = 
0.02. For a value of P c ,3,m above 0.02, we expect to detect image of the slot on the recording medium if 2-2T sensitivity is 
demonstrated. 

Fig. 12 provides contrast S pc3jM for the chosen film (20 lp/mm) for part thickness of 2.4 mm. 
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Fig. 12: Simulated equivalent contrast for film, detector resolution 20 lp/mm, source size 0.2 mm 


Film provides better resolution, but for lower slot depth, the film contrast is low. The film response is also function of 
the film density. Therefore, a different film density response function can be obtained by changing the exposure. A base 
density / 0 of 0-45 was used in this simulation. A higher base density e.g. 2 to 3.5 is preferred for better contrast. 


Fig. 13 provides contrast C pc>3> M for the chosen detector for the part thickness of 2.4 mm. 


X-ray Angle = 0 degree, Source Size = 0.200 mm, Detector Size = 0.125 mm 


Fig. 13: Simulated equivalent normalized contrast for a digital detector with resolution 4 lp/mm 



Slot Depth, mm 


The digital detector provides coarse resolution compared to the film but provides uniform contrast sensitivity compared 
to the film. 


17. USE OF X-RAY PARAMETERS AND SIMULATED CONTRAST 

The X-ray parameter can be used in the POD studies. The crack widths should be measured using a microscope. 
Destructive microscopy on a few cracks will prove to be useful in estimating typical crack widths and crack length-to- 
depth ratio. The X-ray technique sensitivity is assumed to be controlled and maintained at its nominal value using 
measurements on a step wedge and a resolution standard. Experimental correlation of the signal response with the X-ray 
parameter ( P c or P cM ) or pass/fail detection result for a range of X-ray parameters should be obtained. Contrast-to-noise 
ratio should be noted. Image processing can help in improving flaw detectability depending upon the image processing 
routines used. Therefore, a chosen processing routine should be applied to the entire dataset during a POD study. A POD 
curve can be established with the X-ray parameter as the independent variable. The results may establish an a90/95 in 
terms the X-ray parameter for a given, X-ray technique and material. 
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Alternately, the simulated contrast (S) or the normalized contrast (C) can be calculated and correlated with the measured 
contrast ( S and/or C) or with the flaw detection pass/fail result. A POD curve can be established with the simulated 
contrast ( S or C) as the independent variable. 

18. COMPUTER APPLICATION 

A computer application or calculator called “X-Ray Parameter” is developed based on the first 3 models. An example of 
MTF is simulated within the application. It allows choice of detector line-pair resolution. The detector calibration 
sensitivity k can be input along with the source -part-detector geometry. Fig. 14 provides a picture of the main panel of 
the application. The application provides calculated values for the shadow widths, X-ray parameters, MTF values and 
contrast. The data and results for a 70% deep slot with X-ray angle of 5 degree are shown in Fig. 14. 



Fig. 14: X-ray parameter application with example of a 70% deep slot 



Fig. 15: X-ray parameter application with example of a 2-2T slot 

Fig. 15 shows results of a 2% deep slot that is 2% wide similar to a 2-2T penetrameter with X-ray angle of 5 degree. The 
X-ray parameter and contrast values are about same between Fig. 13 and Fig. 14 calculations. The comparison indicates 
that the 70% slot provides a contrast value equivalent to that of a 2-2T penetrameter hole in this case. In the above 
example, the model 2 and 3 outputs are comparable. 
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19. CONCLUSIONS 


This approach assumes that the X-ray parameter and technique requirements including the part material requirements, 
the detector response calibration characteristics together provide a complete description the X-ray flaw detection 
application and provide sufficient information for assessment of the NDE capability. Although, the X-ray scatter is not 
modeled, it is assumed to be low, uniform and controlled by the technique, material, and detector requirements. The 
paper provides derivations of expressions for the X-ray parameters, normalized contrast, contrast, and shadow widths in 
an X-ray set-up for plate inspection. It provides derivations that establish correlation between the X-ray parameters and 
contrast; and the X-ray parameter and normalized contrast. An approach to establish a90/95 for the X-ray parameter is 
provided. A process of obtaining the detector response curve as a function of ratio a/t using a step wedge is described. 
Detector response curve is used in computation of the simulated contrast and simulated normalized contrast. The 
simulated contrast and normalized contrast can be used for establishing correlation with actual contrast and flaw 
detectability. A computer application, given here, provides a simple way to evaluate whether an X-ray technique is 
capable of detecting given size flaws reliably by providing values of the X-ray parameters, contrast and normalized 
contrast that could be compared with corresponding a90/95 values that have been established for similar applications. 
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